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Abstract. Different hydrodynamic models of the heliospheric interface have been presented mean- 
ON \ while, numerically simulating the interaction of the solar wind plasma bubble with the counter- 

streaming partially ionized interstellar medium. In these model approaches the resulting interface 
flows are found by the use of hydrodynamic simulation codes trying to consistently describe the 
dynamic and thermodynamic coupling of the different interacting fluids of protons, H-atoms and 
pick-up ions. Within such approaches, the fluids are generally expected to be correctly described 
by the three lowest velocity moments, i.e., by shifted Maxwellians. We shall show that in these 
approaches the charge-exchange-induced momentum coupling is treated in an unsatisfactory repre- 
sentation valid only at supersonic differential flow speeds. Though this flaw can be removed by an 
improved coupling term, we shall further demonstrate that the assumption of shifted Maxwellians 
in some regions of the interface is insufficiently well fulfilled both for H-atoms and protons. Using a 
Boltzmann-kinetic description of the proton- and H-atom- distribution functions coupled by charge 
exchange processes we emphasize the fact that non-negligible deviations from shifted Maxwellians 
are generated in the interface. This has to be taken into account when interpreting inner heliospheric 
measurements in terms of interstellar parameters. 
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INTRODUCTION 



The mutual interaction of plasma and H-atom gas flows in the heliosheath, in view 
of large Knudsen numbers K n = X ex /L , needs a kinetic treatment of charge-exchange 
induced coupling processes with typical mean free paths X ex larger than typical structure 
scales L, as was already emphasized by Osterbart and Fahr (1992) or Baranov and 
Malama (1993). In kinetic approaches the distribution function of the H-atom gas needs 
to be described by a Boltzmann-Vlasov integro-differential equation (see, e.g., Ripken 
and Fahr, 1983; Osterbart and Fahr, 1992; Baranov and Malama, 1993; Pauls and Zank, 
1996; Fahr, 1996; McNutt et al. 1998,1999; Bzowski et al. 1997, 2000). 

Generally spoken, for time-independent problems one would have to start from the 
following typical Boltzmann equations: 

(v-V r )fi+(F-V^f i = f j (r,v)J /,ft?)v w/ (?,?)a(v«/)dV- 

f 3 
-fi(r,v)J fj(r,v)v re i(v,v')o(v rd )d 3 v (1) 



where indices i, j can be used to denote /#(r,v) and f p (f,v) as the velocity distribu- 
tion functions of the H-atoms and the protons, respectively, r and v are the relevant 
phase-space variables, F are forces per mass, v re i denotes the relative velocity between 
collision partners of velocities v and v', and <r(v re /) is the velocity-dependent charge ex- 
change cross-section. Due to the fairly laborious mathematical tractability of the above 
Boltzmann equation, many authors have preferred to change over from Equ.© to a set 
of hydrodynamic moment equations (for a review see Zank, 1999) thereby only admit- 
ting for the lowest moments of the distribution function to be different from zero, i.e. 
density n ; , bulk velocity £/,•, and scalar pressure niKTi = Pj. For the case of stationary 
ionospheric plasma-gas couplings this had led Banks and Holzer (1968) to the following 
equation of motion: 

Q • V r ) Q + J- (V, ■ *,) -</}) = *L = _^^L (2 ) 

/ ntnii nirrii m t Z>, ; 

where 'Pj denotes the pressure tensor, Cy = Uj — U{ is the differential drift velocity 

between the fluids i and j, (Pi) is the velocity-average of the external forces, Ay are 
forces due to charge-exchange induced momentum transfers, and D ;/ is the ion-neutral 
diffusion coefficient given by: 



One should, however, clearly keep in mind, that the above form of a macroscopic 
equation of motion was derived under the simplifying assumption that Qj/ ' ■sJlKTijmi <?C 



1 (i.e., that highly subsonic differential drifts prevail) and that the charge exchange cross 
section o ex can be taken as independent of velocity. 



SUPERSONIC AND QUASI-SONIC DIFFERENTIAL DRIFTS 



Since in the heliospheric interface, due to the fact that here Qj/ 'y/2KTi/mi > 1 prevails, 
these above made two assumptions are not fulfilled and charge exchange coupling of 
the two fluids "i" and "j" has to be treated in a different manner. Again using shifted 
Maxwellians with isotropic temperatures T p and 7// (as done by Holzer, 1972 ; Fahr, 
1973; Holzer and Leer, 1973; Ripken and Fahr, 1983; Isenberg, 1986; Fahr, 1996; Lee, 
1997) and taking velocity-independent cross sections a ex then suggests to present the 
above mentioned charge exchange momentum coupling term Au in the following form: 

Ay = O re i (v r el) mninjipj - £/;) = TQj = Q (4) 

with (v re i) being the double-Maxwellian average of the relative speed between protons 
and H-atoms as given, e.g., by Holzer (1972) in the form: 



<->=i/f6+l>+(«)~ <« 



The problem appearing with this apporach when treating the passage of neutral in- 
terstellar gas (LISM H-atoms) through the plasma interface ahead of the solar system 
can easily be identified: The problem essentially is comparable to the passage of an 
H-atom gas flow through a predetermined quasistatic plasma structure simulating the 
region downstream of the expected outer interstellar bow shock and ahead of the stag- 
nation point at the heliopause. In a one-dimensional approach for the region along the 
stagnation line (z-axis!) this LISM plasma ahead of the heliopause, due to its very low 
sonic Mach number, can be taken as quasi-incompressible and nearly stagnating. To de- 
scribe the charge exchange imprint of this pre-heliopause plasma sheath on the H-atom 
flow at its penetration through this wall one traditionally uses the following set of equa- 
tions (see Holzer, 1972): 
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where y is the polytropic index taken as identical for both protons and H-atoms. As 
seen from Equ.©, the H-atom mass flow is constant yielding Cq = PqUho = PhUh- In 
addition with introduction of the normalized space coordinate E, defined by z = E,D (D 
being the linear extent of the plasma wall and the quantity A = D/X = Do re in p (E, = 
and E, = 1 mark inner and outer border of the plasma wall) one then obtains the following 
characteristic equation (see Fahr, 2003): 



d_ _ V rel A(A p P p -P H + ±C U H (y+l)) 
dE, H yPh-CqUh 



(9) 



where A p = p p /pH is used. 

Starting the integration at E, = with supersonic H-atom inflow velocities, i.e., with 
U H0 > JPho/Pho > one first obtains physically meaningfull results for Uh and P# with 
increasing values of E,. At a critical point E, = E, c > 0, however, where locally the 
equality jPhc — CqVhc is reached, the integration of the upper system of differential 
equations cannot be continued, since a singularity of an O-type appears which cannot 
be avoided (see, e.g., Kopp and Holzer, 1976). This "neuralgic" point can, however, 
be eliminated when instead of Equ.© a more refined expression for the term of the 
charge-exchange induced momentum exchange between plasma and H-atom flow is 
used which considers the velocity-dependence of o ex = (A+51og(v/vo)) , A and B 
being constants, and the individual relative velocities v re /(v) as, e.g., carried out by 
Williams et al. (1997), McNutt et al. (1998, 1999) or Fahr (2003). As shown by the 
latter author, the corresponding expression, which is valid for moderate and small Mach 
numbers M#, while Equ.© is only justified for large Mach numbers, can be brought 
into the following form: 



Qsub = n 



-giM H - Vk (-9giM H - 2g 2 M H + 5agiM H + 2ag { M H ) 
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(10) 



Here it should be noted that, compared to the above expression, the usually applied 
expression used in Equ.Q, is only justified for the case of high Mach numbers Mh, and 
when written in a manner analogous to Equ. dTO"!) using the above introduced quantities, 
attains the following form: 



Qsuper = -UM H J-^-(l +a) + aMJj (uh/uh) (11) 

In the above expressions the following notations have been used: a = Th/T p ; M h = 
PhUh/yPh; gi = T 5 (l + ^);g2 = gi -V^^;n = ^n p n H o rel ^2KT p /m^/2KTH/m. 
This evidently means that depending on prevailing Mach numbers Mh the plasma-gas 
friction force, based in the past on Holzer's term, was either over- or underestimated in 
published hydrodynamical theories (see Fahr, 2003, for details). 

With the newly derived expression one obtains, instead of Equ.©, the following 
characteristic equation: 

dU H __ Mrei {^pPp-PH - \CqU h (7- 1)) + 7DUhQiM u h,Ph) 
d$ jPh-CqUh 

which now has an X-type critical point with an avoidable singularity condition requiring 
simultaneous vanishing of the numerator and the denominator. Solutions of the set of 
Equs.® through © with application of the newly derived expression (fT0"l) instead of 
(fTTTl are shown by Fahr (2003). 



CONCLUSIONS AND OUTLOOK 

It is clearly manifest from the study of works by Baranov and Malama (1993 ), Baranov 
et al. (1997), Fahr (2000), Fahr et al. (2000), Mueller et al. (2000) or Izmodenov 
(2000, 2001) that the locally prevailing Mach numbers Mh of the relative flows between 
the LISM proton plasma and the LISM H-atom fluid in the heliosheath region (i.e., 
post-bow-shock Mach numbers Mh) generally are found to be smaller than Mh = 2. 
Consequently and strictly speaking, in this Mach number range, the adequate momentum 
exchange term for hydrodynamic approaches must be taken in the newly derived form 
given by Equ. (fTOl . instead in the conventionally used form given by Equ. (fTTTl . Since 
the effective momentum exchange rate described by this new expression, depending on 
prevailing local Mach numbers is smaller or greater than that taken into account by the 
conventionally used term, one can presume that the adaptation of the LISM H-atom flow 
to the bow-shocked LISM proton plasma in the interface region ahead of the heliopause 
operates differently from what is described up to now. 

Hereby the question how to apply the above derived momentum exchange expression 
(fTOl - even though the real interface plasma is not stagnating, but is a two-dimensional 
plasma flow with locally variable properties - is relatively easy to answer: Imagine that 
the plasma flow and the H-atom flow in the interface in nearly all simulation codes con- 
ventionally first are calculated without taking into account the effect of charge exchange 
interactions by only integrating the relevant hydrodynamical differential equations for 



mass, momentum, and energy flow conservation wihtout coupling terms on an appropri- 
ate spatial grid system. Then in a second step of the integrations as usually practiced at 
all grid points the resulting charge exchange interaction terms are evaluated, describing 
local exchanges of momentum and energy per unit of volume and of time. Instead of 
using the conventionally used momentum interaction term given by Equ. lfTTT) . one now 
applies the newly derived term given by Equ. dTOb which when evaluated in the local 
rest frame comoving with the local plasma bulk flow attains just the form identical with 
Equ. (fTOt and thus given by: 

Q* ub = IIM* H (l gl -^(-9g 1 -2g 2 + 5ag l +2ag l M* H 2 )) (13) 



besides the fact that M* H now has to be taken as defined by: M* H = ( U p — Uh ) / \/2AT#/m. 

Since forces are invariant under Galilean transformations (i.e., U p <^ c), one can 

simply, in the next run of integrations, now take into account the ad-hoc forces Q* ub — 

Qsub to remodel the H-atom fluid, and —Q su b t0 remodel the proton fluid. This procedure 
as usual can then be iterated until convergence is achieved. In a future work we are 
going to apply this above mentioned procedure within the frame of our existing five- 
fluid simulation code to better model the interface flows (see Fahr et al., 2000). 

An additional flaw in the hydrodynamic simulation codes is due to the assumption of 
highly relaxated distribution functions f p and fy in the form of shifted Maxwellians. 
This assumption has been shown to be not sufficiently well fulfilled in the heliospheric 
plasma interface since rapid charge-exchange-induced injections of new particles per- 
manently keep the resulting distribution functions away from relaxed hydrodynamical 
ones (see Fahr and Bzowski, 2004). Even the proton distributions in some regions of 
the interface develop pronounced non-equilibrium features, especially in regions where 
the proton densities are low and proton temperatures are high, i.e., where Coulomb re- 
laxation processes have typical periods x pp larger than the injection periods x ex . While 
we refer for details to the paper by Fahr and Bzowski (2004), we here may give helpful 
estimates to characterize the resulting deviations to be expected: Representing the actual 
proton distribution function by f* = fo p + f\ p , with fo p being the hydrodynamic "core 
distribution" delivered by the hydrodynamic multifluid code, one can integrate the total 
rate 8n\ p of protons relaxing by means of pp-Coulomb collisions per unit of time and 
volume towards the core distribution by: 



8n lp ~ -L [\f p -f 0p )d\ = £* (14) 

*PP J pp 

where n\ p is the total density described by f\ p . To achieve stationary conditions, this rate 
must just be balanced by the charge exchange injection rate thus yielding the following 
relation: 

C,P = — - - nQHOex (vrel,p,H) ?pp = — 05) 

n Qp tex 

In the outer interface, where x pp is of the order 2 ■ 10 7 s while x ex is or the order of 
2 • 10 9 s (see Fahr and Bzowski, 2004), one thus obtains q p ~ 10~ 2 . In the inner interface 



inside the heliopause, however, where x ex ~ 4 • 10 8 s and x pp ~ (T p /n p ) ~ 10 12 s 
should be valid, fairly strong perturbations of the relaxated distribution functions must 
be expected. 
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